Squeezing generation and revivals in a cavity-ion system in contact with a reservoir 
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We consider a system consisting of a single two-level ion in a harmonic trap, which is localized 
inside a non-ideal optical cavity at zero temperature and subjected to the action of two external 
lasers. We are able to obtain an analytical solution for the total density operator of the system and 
show that squeezing in the motion of the ion and in the cavity field is generated. We also show that 
complete revivals of the states of the motion of the ion and of the cavity field occur periodically. 



I. INTRODUCTION 

In the last two decades much attention has been given 
to the study of squeezing both theoretically and exper- 
imentally. Squeezed states may be used to improve the 
signal-to-noise ratio in optical communications [1| and in 
spectroscopic experiments 2, 3]. There is hope that it 
can also be used in very sensitive experiments, like the 
detection of gravitational waves jj, iSj] . 

Research on squeezing has been done in several con- 
texts. Here we are particularly interested in systems 
of atoms or ions inside an electromagnetic cavity. It is 
well known that squeezing may be obtained in the in- 
teraction of one mode of the cavity with a two-level or 
a three-level atom. The one photon and the two pho- 
ton micromasers and microlasers are examples of 
such systems. Selective atomic measurements in cavity 
Q.E.D. have been proposed to enhance squeezing (up to 
75%)in the Jaynes-Cummings model 8]. Villas-Boas et 
al. [9j proposed squeezing an arbitrary radiation field 
state previously prepared in a high-Q cavity by the dis- 
persive interaction of a single three-level atom with a 
classical field and a cavity mode. Massoni and Orszag 
[To| introduced a novel way to produce squeezed light in 
an optical cavity through the transferring the squeezing 
from the motion of a three-level trapped ion to the cavity 
mode. Two external laser fields drive the atom and the 
squeezing in the motion is generated by an external elec- 
tric field. Here we have considered the simpler case of a 
two-level ion trapped inside the cavity. This system was 
originally considered by Zeng and Lin and developed 
by many authors ^ei. an analytical solution 

of the master equation was obtained in the case that the 
ion is inside a bad cavity and is excited by only one laser 
field. In this case the energy may be transfered back and 
forth from the cavity field to the motion of the ion. 

In this paper we discuss a scheme for generating 
squeezing both in the motion of an ion, driven by two ex- 
ternal laser fields, and in a mode of the electromagnetic 
field inside a non-ideal cavity at zero temperature. Even 
though we consider a non-ideal cavity, we are able to ob- 
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tain an analytical solution for the total density operator 
at any time and show that the motion and the cavity field 
are always squeezed in one of their quadratures if the ratio 
between the intensity of the lower and higher frequency 
fields is less than one. Remarkably this system presents 
revivals, the vibrational state periodically turns into the 
same squeezed vacuum state and the cavity field periodi- 
cally visits the vacuum state. Although we are concerned 
here to a specific system, our analytical results may be 
useful in any situations were the effective Hamiltonian 
corresponds to a general linear coupling of two oscilla- 
tors, one of them being connected to a reservoir at zero 
temperature. In Sec. II we present the model that couples 
the ion to the cavity field. In Sec. Ill we give an analyt- 
ical solution for the master equation showing explicitly 
the revivals and study the generation of squeezing. In 
Sec. IV we summarize our results and present our con- 
clusions. In the Appendix we present a detail calculation 
for obtaining the solution presented in Sec. III. 

II. THE MODEL 

The system we are considering is very similar to the 
one discussed in Ref. A two-level ion is bounded 

by a linear Paul trap, which is inside an optical cavity. 
The ion has mass m and the effective potential created 
by the Paul trap may be approximated by {1 / 2)mv'^ SP' , 
where x is the displacement operator of the ion from its 
equilibrium position and v is the classical frequency of 
oscillation of a particle in this potential. The trapping 
direction x coincides with the axis of the cavity. The 
electronic levels of the ion, |e) and \g), are separated by 
the energy Tiljq and are quasi resonant to a stationary 
mode of the cavity of frequency lJc- This system is then 
described by the Hamiltonian H : 

H = Ho + Hi , 
Hq ~ hujQ\e){e\ -\- huJcO-''a 
Hi — hgcShi{kcx)a\e){g\-{']i.c. , (1) 

where a and a) (6 and w ) are the annihilation and cre- 
ation operators of the cavity mode (vibrational quanta) . 
In Hi , fcc is the cavity wavevector and Qc is the ion-cavity 
coupling constant. We have taken a sinus function as the 
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FIG. 1: Level scheme of the ion. The excited level |e) and the 
ground level \g) are coupled by laser fields of frequencies lji 
and UJ2 and one cavity mode of frequency cOc, detuned by A. ojq 
and u are the atomic transition and the vibration frequencies, 
respectively. 



ion starts in its ground state, we may substitute the op- 
erator (Tz by — 1. In this case, the effective Hamiltonian 
Q may be rewritten as 

Fcff = ihfli {a^b - aP) + inQ.2 (a^S^ - ab) . (4) 

The first term of Eq. describes processes where 
the sum of the number of photons and of vibrational 
quanta are conserved while the second term describes 
processes where the same number of photons and vi- 
brational quanta are created or absorbed simultaneously. 
These processes allow the exchange of quantum informa- 
tion between the cavity field and the vibrational motion 
and the generation of squeezing. 

III. THE DENSITY OPERATOR 



cavity standing wave mode and set the minimum of the 
trapping potential at a cavity node. The eigenstates of 
Hq are tensor products of the electronic states, |e) and 
times Fock states \n)c and |to)«, associated to the 
cavity field and vibrational quanta. The position opera- 
tor is related to the operators b and ¥ hy x ~ 5x{b + ¥), 
with 5x = yjh/(2mv) being the uncertainty of position 
in the vibrational ground state. 

Now we let two external lasers act on the ions. The 
laser frequencies, lui — lOc — v and uj2 — lo^ + v, and 
the detuning, A = luq ~ lOc, are chosen in order that 
Raman transitions among the level \g)\n)c\m)v and the 
levels \g)\n + l)c|''Ti ± l)v, Iff)!"- — ^)c\'m ± l)t, may occur 
(see Fig. ??). The total interaction Hamiltonian that de- 
scribes the coupling of the internal and external degrees 
of the ion with the cavity and with the laser fields, in the 
interaction picture with respect to Hq, may be written 
as 



JAt 



Hint{t) = h{gcsm{kcx{t))ae'^" - igie 



(2) 



where fci , k2 are the x components of the laser wavevec- 
tors and 51,(72 are the ion-laser coupling constants. x{t) 
is the position operator in the interaction picture. 

We consider that the Lamb-Dicke limit {rjc — kcSx <C 
1, 771 = ki6x ^ I, 772 = k2Sx ^ I) is vahd. We 
also assume the rotating wave approximation and that 
Wo > A > r]igi,r]2g2,ricgc,v. For t > 1/A, we obtain, 
following the usual procedure for the adiabatic elimina- 
tion of one of the electronic levels^^, the effective inter- 
action Hamiltonian: 

V^'g = —ihfli (a'^b — aP)a;, — ihn2 (a^6^ — ab)(Tz 
+ higl+gl)/Aa,, (3) 

where ct^ = |e)(e| - \g){g\, fli = ??cffiffc/A and 
^^2 = 'rjcg2gc/ A- The third term is the usual Stark shift 
and will be incorporated in Hq by redefining wq- Non 
linear terms in rj will be neglected by assuming that the 
average number of excitations are not very large. If the 



Sources of dissipation in the mechanical movement in- 
side a Paul trap may be attributed to chaotic fields due to 
strain fields in the trap electrodes and to the electronic 
decay of the upper level |e) 0, 0|. They are much 
easier to control than the losses of the cavity field. For 
example, heating in recent Paul tra p e xperiments have 
been limited to 1 quantum per 6ms Also the heat- 
ing due to the electronic decay could be made small if 
the detuning is sufficiently large. Here we only consider 
cavity losses and neglect the interaction of the ion vibra- 
tional motion with the environment. Assuming a linear 
interaction of the cavity field with a reservoir at zero tem- 
perature and using the Born-Markov approximation we 
may write a master equation were the losses are included 
via a Lindblad form [13 



^ ^[%srPT{t)] + cpT{t) , 

where V^ff is the effective Hamiltonian (0} and 



(5) 



LpT — ■^(SapTfl^ — a^apT — pra^a) . (6) 

In the Appendix, we solve this master equa- 
tion when the initial state is the vacuum state, 
prp{0) = |0)c c(0| |0)^, v{0\ . We obtain for the density op- 
erator at time t 

00 oc 
m— n— 

X Q"rinc{t)Ac{t))Q"rinv{t),Ut)), (7) 
where the functions f(t) and g{t) are given by 



fit) 



cos(At) + ^ sin(At) ] e-"^*/" 



g{t) = :^sin(At)e-^*/\ 



(8) 



with A^ 
with a = 



- nj-nl- 7VI6. The operators 
c,v, labeling the cavity field and vibrational 
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motion respectively, are defined as 



= E '^TiO SAO i?™+"-'='^ri) ^.(0^ , (9) 



fc=0 



with S'o-(^) being the squeezing operators that act on the 
cavity or in the vibrational states: 

The coeiEcients C™'"(^) are given by 



(10) 



C7A0 = 

min(n,fc) 

X E 



^m,r! /^^ _ /(m + n - A;)!fc! 



mini 
ml 



{k-l)l{m-k + l)ll\{n-l)l 



max(0,fc— m) 

X (coshO""''+^'(sinhO"+''"^' 



(11) 



and the operators i?™'", which satisfy the relation 
_R"^™ = (_R™'")^, are defined for m > n as 



ilkl 



X P 



fc,fc- 



fc=0 



n + 1 



i!(fc + m-n)! (n + 



I fc + m — n) 



where 

k 



(j-0!(fc-j)! j! ' 



(12) 



(13) 



J— max(0,/) 

may be identified as a Jacoby polynomial. The functions 
nc(t), nv{t), S,c{t) and £,v{t) in Eq. are given by 



n,(t) = -1/2 + v/1/4 + i^At) - (?' - 1) ''AtV , 

If ^/M 1 1 /^l/2+_^i+lK(t)\ 

^ i^"l l/2-(g-l)..(t) j' " = 

Cc(i) - -iec(t)i , ut)^\ut)\, (14) 



where q = fli/il2, i^c{t) — g{t)'^ and 

0^ — 1 



(15) 



The reduced density operators are obtained by taking 
the partial traces of Eq. ((T)). By using the property (see 
Eq. lU26t l. 



tr,Q™'"(n,e) =^,„,o^„,o, (16) 

we easily get 

Mt)=Q''.°inAt),Ut)) , a = c,v. (17) 

In the Appendix we show that Q^^{n, ^) has the form of 
a "squeezed thermal state" (see Eq. (jA41|) '): 

oc 

Q?(n, = E (f, + i)k+i SAO \k)aAk\ SACV ■ (18) 



/c=0 



A. > fil : squeezing generation and revivals 

When > ri^, the system reaches a steady state. 
From Eqs. (P) . itTH) . itTsll . we see that in the limit t 
oo, ny{t),nc{t),^c{t) 0. Therefore the total density 
operator in the steady state may be written as 



with 



where 



Pt =|0)cc(0|/5r, 
pf = 5„(C) |0).„(0| S,0, 



-1 / r>i + »2 
^ 2\\ni-n2\ 



(19) 



(20) 



(21) 



Eq. l(T^ shows that when the system reaches the steady 
state the energy transferred from the lasers to the cav- 
ity dissipates away while the vibrational state becomes a 
squeezed vacuum state. 

When fll > 7^/16 the function f{t) vanishes at 

times 



1 f -7/4 

T„ = — arccos 



"A' 



(22) 



with 7r/2 < arccosx < tt and n ~ 0,1,2,.... Since 
fi'Tri) — Oj the cavity field and the ion motion are in- 
stantly decorrelated at i = t„ : 



/5t(t„) = Pc{T„)pv{Tn) , 



(23) 



as can be seen from Eq. ((TJ. Also at these times ni,(T„) = 
and ^^,(t„) = From Eq. ^ we obtain 

Pv{Tn) = S,,{0 |0)„ „(0| 5„(0t = , (24) 

that is, the reduced operator that describes the vibra- 
tional system returns to the same state periodically and 
this state coincides with the reduced operator in the 
steady state. Notice that for these finite times the cav- 
ity state remains a "squeezed thermal state" with fic(T„) 
varying as r„ increases. 

There is also recurrences in the state of the cavity field 
at times 



TT 

"a 



n = 0,1,2,. 



(25) 



when ^ci^n) — ^c{Tn) — 0, since g(T^) — 0. Therefore the 
cavity field and the motion arc also instantly decorrelated 



at times r^, 



PT«) = /5c«)/5„«) , 



(26) 



and the cavity field returns, with a period tt/A, to the 
vacuum state, 



PcK) = |o). 



(27) 
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For small values of ^^2/^1, the "squeezed thermal state" 
given in Eq. ifTHl) . is very close to an ideal squeezed state. 
In fact, the maximum value for fic(t) and ny{t) is given 
by 



v/l-(r!2M)2 



- 1 



(28) 



which is approximately 0.25(il2/rii)^, for (02/f2i)^ <C 1. 

Having obtained the solution in the case that the initial 
state is the vacuum, Eq. ((TJ, we may easily find analytical 
solutions of the master equation when the initial state is 
the product of coherent states: 

^^■'"(0) = D,{a)\0), MDcia)^ D,{m)v .(0|i).(/3)^ , 

(29) 

where 



(30) 



In this case we have 



(31) 

where px (t) is the density operator for the initial vacuum 
state, which is given in Eq. lO, and u{t) and v{t) are 
defined as 

u{t) = ah{t) + {(3ni/n2 + (3*)9it), 

vit) = {~ani/n2 + a*)git) + (3 fit) . (32) 

The functions f{t) and g{t) were already given in Eq. 
and h{t) is defined as 



h{t) = ( cosAi- ^sinAi 



-7t/4 



(33) 



That is, the solution of the master equation, when the 
initial state is the vacuum displaced by a and P in the 
phase space of the whole system, corresponds to a dis- 
placement of u{t) and v{t) of the solution of the master 
equation when the initial state is the vacuum. Then the 
reduced density operators may be written as 



(t) = bc{u{t)) p,{t) b,{u{t))\ 



(34) 



where Pc{t) and Pv{t) are given in Eq. (|17|l . When fi^ > 
rt2, for any values of a and (3, the cavity field and the 
motion of the ion are decorrelated instantly at times r„ 
and r,'j, given in Eqs. H22() and (|25|l . For a — recurrences 
will also occur at these times, since when i = r„, /(t„) = 
and when t = t'^, 9{Tn) = 0- As a consequence f (t„), 
nvijn), u{t^), fic«) vanishes and^„(r„) = |, Cc«) = 0. 
Therefore these revivals occur, independently of the value 
of f3, at t = Tn to the same squeezed vacuum vibrational 
state obtained before and given by Eq. (|24|l , and at t = 
to the cavity vacuum, Eq. H27|) . In fact, we have verified 



that this last result is valid for any initial state of the 
motion as long as the initial state of the cavity field is 
the vacuum. 

Using Eqs. l(T7|l . ifT^ . ||23J, we may show that the un- 
certainties associated to the quadratures 



a + a' 



Pr 



6 + St 

Y _ p _ 



a — a' 

ITT' 



V2 



zV2 ' 



(35) 



have the form 



^X^tf = (fv(t)-f l/2)e-2«.(*)^ 

AP„{tf - (n,(t) + l/2)e2«"(*), (7 = c,t,,(36) 

in the case that the initial state is a coherent state. Notice 
that AX^(t)2AP^(i)2 = {fia{t) + 1/2)2. 

By substituting the expressions oina{t) and ^ait) into 
the above equations, it is easy to show that: 



AXe(<)2 ^ 


1 


Ql + il2 




2^ 


^2 




1 


- n2 


■9{t)\ 


2 " 


^2 


AX^ity = 


1 


^2 




2 " 


+SI2 


AP„(t)2 ^ 


1 


^2 




2^ 


Qi - n2 



that is, when fJi > 1^2, the quadrature Xv is always 
squeezed for t > 0, while the quadrature Pc is equal to 
1/2 when t = t,j = nir / K and is squeezed otherwise, since 
(1 - f{tf) > for i > and g{tf > for t ^ < . 

In Fig. 2 we show APc{tf and AX.^{tf for Q.2/VI1 
equals to 0.6 and 0.9 and for 7/fii equals to 0.4 and 1. 
The maximum values of the squeezing for the X^ quadra- 
ture are reached at times r„ when the vibrational state 
is equal to the squeezed vacuum state. At these times 
nv{Tn) — and the uncertainty in the quadrature Xy is 
given by 

AX.(..f4e-«- = i2i_5i. ,3a, 

and is equal to 0.25 and 0.05 when Q.2/^1 equals to 0.6 
and 0.9, respectively. The quasi periodicity shown in 
the values of the uncertainties of the quadratures is a 
reflection of the fact that f{t) and g{t) are periodic func- 
tions attenuated by the exponential e~^*/^. The values 
for 7/ili in Fig. ?? are not far from those used in recent 
experiments |23|. For example, if we take .gc/27r = 
6 - 7MHz, 7 = (0.02 - 0.05)gc and 5i/27r « lOMHz, 
A « 5^1 and 77c = 0.2 we have 7/rii ~ 0.4 — 1. 

When < Vl\ <^l+ 7^/16, A = i|A| is imaginary. 
In this case, as t — > 00, the functions cos(A<) and sin(At) 
diverge. However the functions f{t) and g{t), given in 
Eq. ijSJ), do converge to zero and the system still reaches 
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■In 



(VLx + - ^2 cos(2Aoi)) 



{Qi - fi2)(f^i + cos(2AoO) 
uo{t) — a cos(Aot) + a sin(Aoi) , 
vo{t) = /3cos(Aot) +/3sin(Aot) . 



(40) 



where Ag — fl^ 



r^l, and 



/3f^i)/Ao , 
-ar2i)/Ao. 



(41) 



When Aq > 0, \ip{t)) is a periodic function of t with 
a period 27r/Ao. When = rmr/ (2Ao), m = 0, 1, 2, . . . , 
no{tm) = and the state of the motion and the state of 
the cavity field disentangle and are given by: 

l^i) = D,{a)\0)MP)\0),, 

IV-s) = 4(-a)|0)A(-/3)|0)„, 

|V4> = D,i-a) Sc{-0 |0)c SviO |0).„ , (42) 

where f is given by Eq. (|21|l . 



B. < f7i 



FIG. 2; APc(t)^ and AX„(f)2 as a function of the adimen- 
sional time flit for ^2/^1 equals to 0.6 and 0.9 and for ^/^i 
equals to a) 0.4 and b) 1. The squeezing e~^^ is equal to 0.25 
and 0.05 for n2/f2i equals to 0.6 and 0.9, respectively. 



the steady state given in Eq. p9|) . although it does not 
present revivals. 

It is expected that recurrences of the whole system will 
occur when dissipation is absent. In fact, one should be 
able to diagonalize the quadratic Hamiltonian given by 
Eq. Q through a Bogoliubov transformation and show 
that this system is equivalent to two independent har- 
monic oscillators and therefore revivals will occur for any 
initial condition. The remarkable result that we have 
shown above is that recurrences in each individual sys- 
tem occur when dissipation is present. For comparison 
we give below the time evolution of the system when dis- 
sipation is absent and the initial state is the product of 
the two coherent states |V'(0)) = Dc{a)\0)cDy{f3)\0)v : 



(39) 



where 



noit) = i (^Y/l + (r!2/A2)sin2(2Aot) - , 



When ill — the analytical results obtained in 
Eas. H5ll8l) are still valid, but the total system does not 
reach a steady state since it gains more energy than it is 
able to dissipate. When < fij, Cc(0 ^^'^ ^v{t) con- 
verge, as t ^ 00, to the finite value given in Eq. H21|l . 
However nc{t), nv{t) and f{t)g{t) increase monotonically, 
so that the ion-cavity system does not reach a steady 
state and remains always correlated. All uncertainties 
increase monotonically with 172/^^1- 

When ill = fl2, the system still remains correlated all 
the time, since f{t)g{t) never vanishes. Also, both, the 
cavity field and the ion motion, do not present either 
revival or squeezing. The uncertainties in the quadra- 
tures APc(i)^ and AX„(t)^ are equal to 1/2 during all 
the time, as we can see from Eq. (|37ll . while the other 
uncertainties are given by 

2 7^ 
1 1 fi02 

AP^tf = L + ^i^t/2 + e-^'/'-lf. (43) 

Notice that Eq. H21(l is no longer valid for il2/i^i = 1, as 
we may not interchange the limits Q2 ^ ^1 and t 00 
in Eq ifT^ . 

In Fig. ?? we show the uncertainty AX„(t)^ as a func- 
tion of ri2/rii for several times and for a fixed value of 
7 — OAili. Notice that all curves cross at ^1.2/^1 = 1, 
when AXy{t)'^ = 1/2, and increase very fast for ^2/^^! > 
1. 
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FIG. 3: AXv{t)'^ as a function of ^^2/^1 for flit equals to 1 
(dotted line), 5 (dashed line) and 10 (solid line). y/Q,! — 0.4. 

IV. CONCLUSIONS 

We were able to provide an analytical solution for the 
total density operator of the system of a two-level ion 
trapped inside a quasi-resonant cavity subjected to the 
action of two laser fields with frequencies lOc^i', when the 
initial state is a general coherent state. The two external 
lasers generate an effective interaction that corresponds 
to: a) a parametric excitation of the vibrational mode 
and the cavity field (term proportional to fl2) and b) an 
exchange of energy between them (term proportional to 
ill). When rij > the cavity mode and the ion vibra- 
tional motion are both always squeezed and the reduced 
density operator of each mode corresponds to a squeezed 
thermal state. The maximum value of the squeezing in 
the vibrational mode occurs with a period tt/A. When the 
initial state of the cavity field is the vacuum, the reduced 
density operator of the vibrational motion has "revivals" 
of the squeezed vacuum state given in Eq. (|20|l . The re- 
vivals in the motion state occur when the uncertainty in 
one of its quadratures reaches the minimum value given. 
The cavity field returns to the initial state periodically. 
At times when revivals occur the cavity and the motion 
disentangle. 

When fli < ilj the system does not reach a steady 
state and does not present either revivals or squeezing. 

The analytical solution for the total density operator 
may be useful in problems involving the interaction be- 
tween two harmonic oscillators in the case that one of 
them is connected to a reservoir. 
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APPENDIX A: SOLUTION OF THE MASTER 
EQUATION 

In Refs. 0, 1^ the authors have obtained a solution 
for the master equation ((SJ in the case = 0, by using 
ladder operators and an expansion in the eigenstates of 
the Liouvillian superoperator. We were not able to find 
a simple similar solution in the case Q2 7^ using that 
method. Here we present the solution of Eq. jsj using a 
different approach. 

We are interested in finding the solution of the mas- 
ter equation when the initial state is the vacuum, 
Pt{0) — |0)cc(0||0)t,^(0|. It may be written formally as 
follows: 

PtW = e'C*|0)ec(0||0)„,(0|, (Al) 
where K. is the Liouvillian superoperator 

IC = ^(Kff. - -Kff) + J(2a..at - ata. - .a^a) , (A2) 
in 2 

with Vcs being the effective Hamiltonian given in Eq. Q). 

We have adopted the notation A. {.A) for superopera- 
tors that represent the simple action of an operator A to 
the left (right) on the target operator, p, i.e., A.p := Ap 
{.Ap := pA). 

It will be convenient to define the following superoper- 
ators 



M% 






.a — a 




M'L = a. 


Ml - 


.at. 


Ml 




Nl = 


.b-b. 




M'L = b. 


, Nl - 





which obey the commutation relations 

[M'L,M%] = 1 , [N1MV\ = 1, 
[M'L^Ml] - 1 , [Ml, Ml] = 1 , (A4) 

while the remaining relations are null. The Liouvillian tC 
may then be written as 

IC^IC' + n2{MlMl+MlMl), (A5) 

where 

JC' = ~-i/2{MlM'L+MlMl) 
+ niiMlM"! - M1MX) 
+ VLi{MlMl - MlMl) 
+ n2{MlMl + MlMl) 
+ n2{MlM''_ + M''_Ml) . (A6) 

From Eq. ljA3|l it is easy to see that 

X-|0),.(0| = , AAf|0),,(0| =0, (A7) 
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where a = c, u. Thus the superoperator KJ annihilates 
the vacuum: 



/C'|0),,(0||0)„„(0| =0. 



(A8) 



It is possible to relate K, and K! through the transforma- 
tion 



(A9) 



where 
U = 



^\ 1 ^2 

(AlO) 

By exponentiation of Eq. ljA9p we may rewrite the time 
evolution superoperator e'^* as 



(All) 



Therefore the density operator at time t may be written 
as 

p^(fy = e'^e'='*e-"e-'C'*|0)ec(0||0), „(0| , (A12) 

since * does not affect the vacuum state. The super- 
operator *e~^e~^ * may be written as 



(A13) 



To calculate the term in brackets in the above equation, 
we will need the following commutation relations: 



-7/2 







^2 





-7/2 
















Vt2 











Ml 



(A14) 



From Eq. (|A14p and using the Baker-Haulsdorff formula, 
we get 



Ml 



I 



h{t) 






h{t) 

git) 



V 



m Trait) 



^2 



VL2 
git) 

fit) 




ait) git) 



\ 





fit) 



(M%- 
Ml 
Ml 

\Ml 



(A15) 



where the functions fit), g{t) and h{t) are given by 

7 



fit) 

git) 



cos At + 4- sin At e^'^*/". 
4A 



^sinAte"^*/4. 
A 



h{t) = (^cosAt- ^sinAt^e-''*/\ (A16) 

with A^ = nl-VLl--f'^ The 4 X 4 matrix in Eq. 
was obtained by exponentiation of the 4x4 matrix in 
Eq. IIA14p . what can be easily done by noticing that the 
latter may be rewritten in the form — 7/4I-I-AA, where 1 
is the 4x4 identity and A is a matrix such that = — 1. 

Using Eq. (|IlO|l . Eq . lfn5)l and that [U.e'^^Ue'^*^] = 
0, we may rewrite Eq. ljA12|l 



hit) 



- ,>V(t)| 



(A17) 



where 



Wit) = U - e^'^Ue-'^'* 



= VMt), u,{t)) + V,(Ai,(t), v^it)) 

+ Cit)iMlMl + MlMl) , (A18) 



with 



Vai^Ji,v)= ^i{Ml^ ^Ml^)/2^vMlMl. (A19) 

The functions C('^)i Mc(i)i ^^cit), l^-vit), and Vvit), in the 
above equation are given by 



at) 

f-^cit) 
Vcit) 

fly{t) 

Vvit) 



fit)git), 
git)\ 



(A20) 



All the superoperators appearing in the definition of 
yV(t) commute. Therefore we may write 

gW(t) ^ exp[(:{t)MlMl] exp[Cit)MlMl)] 
X exp[VcilJ-cit),i^cit))] 
X exp[VyitJ.yit),iy^{t))]. (A21) 

By substituting Eq. HA21|I into Eq. (|A17p . we obtain 

00 00 

Mt) = EEcw-^" 

m— n— 

X Q'-^-'^ipcit), Mt)) QZ'^iPvit), J^„(t))(A22) 

where 



(A23) 
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with 

0^(^,1.) = e^-(''^'')|0)..(0|. (A24) 
Due the cychc property of the trace we have the identities 
tr,(A^;p)-0 , tr,(AA^p) = , Vp. (A25) 
Then we easily may conclude that 

tr,Q^'"(M,i^) = W„,o. (A26) 

Below we give a useful expression for the operators 
Q'^^{li, v), which will allows us to write the reduced den- 
sity operator Pa-{t) as a product of squeezing operators 
times an operator that describes a Planck's thermal dis- 
tribution. In order to obtain this result we will start by 
defining three new superoperators: 



{Ml ^ +Af^ 2)/2 - MlAfl - M'LAfl 



J^(y — 

which form the closed algebra: 



+ ' 



(A27) 



[B^Xa] - 2X + 2C 



(A28) 



The superoperator Va{^J■^ v) may be expressed in terms of 
these superoperators and, due to the fact that they form 
a Lie algebra, we may decompose e^"'^'^''^^ into a product 
of simpler exponentials. 



^ ^(,A„ ^nBa ^iCa 



(A29) 



where 



= -1/2 +v/(;^+ 1/2)2-^^2 



1 (v^^l^-ji^ 
^ 4 '\z/+l/2 + // 



(A30) 



From Eq. ^\ and Eq. (|X27ll we see immediately that 

e«''"|0)..(0| = |0),.(0|. (A31) 
Also the operator 



i?r(") 



|0>..(0| 



(A32) 



may be calculated by noticing that 



and that 



Then R°'°(n) satisfies 



(X^-7W_^)i?°'°(n) = 0, 

{AfZ - nMl)R°-°in) = . (ASS) 



Using the identities (|A25|I and the definition (|A32|I we 
easily get 



tr,i?:;'"(n) = 1 . 



(A36) 



The solution of Eq. ljA35p under the normalization con- 
dition given in Eq. lfi36ll is a thermal state: 



oo 



lO,0 



(A37) 



Thus the operator Q^ '{n,^) may be rewritten as 

Qy(n,0=e^^-R''f(n). (A38) 
It is possible to sec that 

e«-^" = 5,(0•■^.(0^ (A39) 
where S^iS,) are the squeezing operators 

S.iO = e«/2(^'-^''). (A40) 
By using Eq. l|X37|) . Eq. (|X38l) and Eq. l|X39|) . we obtain 

OO _^ 

= 127^TWTT^^(^)\'')--(''\^-(0^ ^ (A41) 
fc=0 ^ ' 

which is the same as Eq. (|18|) of Sec. III. 

By using Eq. (|A38|) and inserting identities 1 = 

we obtain 



into the definition of Q™'"(n,0, Eq. jJ23)l . 



Q^>"(n,0 = e' 



(e -M+e^^o^^.^^ 



(A42) 

By substituting the equalities 

g-M._A/f^eM. = cosh^TW^ +sinhC7V_^ , 
g-M„_^^gM. = sinh^TW^+cosh^A^J , (A43) 

into Eq. HA42|I we obtain 

= E c';;'"(0 5.(o^"+"-'''N^.(o^(A44) 

fc=0 
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where we have defined 



and 



ATcr n KAU m 



min(n,fc) 



(A45) 



(m + n — ky.kl 
mini 



E 



(k-l)l{m-k + l)ll\in-l)l 

/— max(u,/c— m) 

X (coshO"'"''+^'(smhe)"+''"^' . (A46) 



Using Eq. IIA37|) and Eq. (|A3|) . we obtain the expression 
of for TO > n, in the Fock basis 



ilkl 



1 



fc=0 



TOJ(fc + TO-n)! + 



X P 



k.k — n 



n + l 



|fc + TO - n)c, a{k\ , 



where 



/ iv,-i ij + m)l 

TTTTI 7TT-7 



J— max(0,t) 



(A47) 



(A48) 



The expression of i?™'" for m < n can be easily obtained 
from the property i?™'" = 
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